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1. INTRODUCTION 

Creep has been a topic of investigation for many decades. Meaningful 
results have been obtained by investigators both from the micromechanistic 
approach and the macroanalytical or continuum mechanics approach. A rich 
literature is available and references may be found in the books by 
Garofalo [1], Odquvist [2], and Rabotnov [3], among others. 

Both the micromechanistic and continuum approaches can lead to fruitful 
results and each can benefit from the other. "The micromechanistic approach 
provides knowledge of the processes that control creep and provides guide- 
lines for defining more clearly material properties and for designing better 
materials for specific applications. The macroanalytical approach can pro- 
vide basic relations that are broad in scope and can lead to improved pro- 

ft 

cedures for designing structures." 

The purpose of this report is to discuss further the behavior of creep 
in metals using the approach of continuum mechanics. Several theories have 
been proposed in the literature to describe creep using this approach. How- 
ever, there still exists an unresolved problem related to the role played by 
strain-hardening during creep. This problem is important in the investiga- 
tion of creep subject to variable stress. 

In most investigations, the subjects of stress-strain relation and 
creep are treated separately, so that they do not bear direct relationship 
between them. In the writers' view, the aforementioned fact contributes 
greatly to the difficulty associated with the investigation of creep under 
variable stress. Fortunately, recent progress in the understanding of the 
strain-rate and strain-rate history effect in the theory of constitutive 
* 

The words of Garofalo [1] are being quoted. 
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equation has made it poeelble to deviae a unified approach which would bring 
diversified material behaviors, such as constant-strain-rate stress-strain 
relation, creep, and stress relaxation into a common workframe. 

The present report is written with such thought in mind. This viewpoint 
is shared by such research workers as Valanis and Lalwani [4] and Cernocky 
and Krempl [5]. In this unified approach, creep is viewed as a special case 
of the general mechanical behavior of material. The constitutive creep 
equation is reduced directly from the general constitutive equation of the 
material under study but with the condition of constant stress imposed. 

Since the constant-strain-rate stress-strain relation is also reduced from 
the same general constitutive equation, it is evident that a correlation 
between the two areas of interest can be established. Thus, material con- 
stants and functions determined from the constant-strain-rate stress-strain 
curves will appear without alteration in the creep equation. The number of 
unknown parameters in the creep equation is thus greatly reduced. 

In this report, Valanis' endochronic theory of viscoplasticity [6,7,8] 
is applied to tackle the problems associated with creep. The endochronic 
theory has been previously applied to investigate the creep behavior of 
metallic materials. Valanis and Lalwani [4] developed a nonlinear evolution 
equation for the thermodynamic internal state variables using the concept of 
absolute-reaction-rates theory. The equations are then applied to predict 
creep at moderately large strains under constant stress from the experimental 
data obtained in stress relaxation. Wei [9] investigated creep using the 
Gibbs free energy formulation and the stress-defined intrinsic time. 

The approach taken in this investigation is different from those men- 
tioned above, although all are within the framework of the endochronic 



theory. In particular, the concept of intrinsic time introduced by Valanis 
[8] is employed in this investigation. 

The problems of creep recovery ard stress relaxation are also treated 
using equations derived from the same general constitutive equation by 
imposing appropriate constraint for each case. It is shown in this report 
that the theory agrees quite well with experimental results of Uang and 
Onat [10] for Aluminum 1100-0 at 150°C (300°F). 
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2. BRIEF SUMMARY OF ENDOCHRONIC THEORY 
OF VISCOPLASTICITY 


The endochronic theory of viscoplasticity developed by Valanls [6,7] 
is based on the notion of intrinsic time and the thermodynamic theory of 
internal variables. Since most of the materials are, in general, strain 
history dependent, a time measure dc is defined such that 

d ' 2 ' p ijkt dE ij d£ k« + 8 dt2 (1) 

where is the strain tensor and is generally a function of and 

a positive definite material tensor; g is a material function of c^; and 
t is the real time. In addition, a time scale z(C) is introduced such that 
dz/dC > 0. 

This concept together with the thermodynamic theory of the Internal 
variables gives the following explicit constitutive equation for isotropic 
materials under small isothermal deformation: 


f 


A(z - z') 


3e 


kk 

3z' 


3a, 


dz* + 2 p(z - z’) “ dz' 

Jn az 


(3) 


where a., is the stress tensor, e,, is the deviatoric part of e,,; 6,, is 
ij i j ij ij 

the Kronecker's delta; and X(z) and y( z) are heredity functions. But the 
definition of intrinsic time in equation (1) has led to difficulties in 
cases where the history of deformation involves unloading. Valanis [8] 
has since introduced a new concept of intrinsic time to overcome these 
difficulties. In the one-dimensional case the new Intrinsic time c is 


defined as 



5 


t 


<1C - 


de 



(4) 


where k^ is a positive scalar such that 0 < k^ < 1 and Eg is the elastic 
modulus. Generalizing to three dimensions and r internal variables, a 
strain like tensor 0^ is defined as 

9 ij " e ij " *ijki °ki 

and 



'ij 



( 6 ) 


where is a positive definite symmetric fourth-order tensor; Q. . is 

ijK* lj 

the deviatoric part of 6^; is the shear modulus; and s^ is the devia- 
toric part of o^. Based on the formulation in reference [6], the response 
of metals in the small deformation region with an elastic hydrostatic re- 
sponse can be written as 


kk 


3K e 


kk 


and 


(7) 


8 


ij 



u(z - z') 



( 8 ) 


where K is the bulk modulus. If the function ii(z) is defined as 
y(z) - Wq G(z) 


( 9 ) 
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where G(0) ■ 1, end using Che Leplace Transformation technique, equation (8) 
together with equations (6) and (9) reduce to 


f* dQ ±1 

" 2W 0 J 0 P(Z “ * f) d *’ (10) 

where p(s) is related to G(z) by the integral equation 

p(z) “ k i J P(* " **) dz* « G(z) (11) 

The solution of this integral equation in the case of two Internal variables 
is obtained by the Laplace transform technique under the assumption of 
approaching to 1. The result is 


E q p(z) 


— 5(z) + E. e”“ z 
a 1 


( 12 ) 


where a^, a, and 
In the case 


E^ are material parameters and 6(z) is the delta function, 
of uniaxial stress, the constitutive equation is given by 


-of 

where 0 * e - ~ is 
E 0 

dt » |de| 


p(z - z') dz* (13) 

the plastic strain, and equation (4) can be written as 

(14) 


Because of the complexity of the constitutive equation and the fact that 
the number of material parameters Increases with the number of internal 
variables, it is natural to strive for the minimum number of internal 
variables which will suffice in the accurate and acceptable prediction 
of the material reaponse. Two internal variables are used throughout 
this work. 
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The new form of the endochronic theory has recently been extended by 
Wu and Yip [11] to describe the strain rate and strain rate history effects 
of material behavior. In this connection, it was shown that the intrinsic 
tine is governed by the relation 

dc-k(0)|de| (15) 


where 8 is the plastic strain rate. 

By combining equations (12), (13), and (15) it was shown that yield 
stress a can exist and is related to the strain-rate by 

y 


o (fi) 
y k(e) V*' 


where o y is the initial yield stress at reference strain-rate. Making use 
of the expression (Ref. [6]) 


z - ± ln(l + Be) 


where B is a material constant, the constant-plastic strain-rate stress-strain 
relation was obtained 


o - 0 y (l + B x 9) + (o q - o y )(l + 0 x e>[l - (1 + 3^6) ] 


in which 




n - |+1 


and 0 q is the Intercept of the asymptotic line of equation (18) with the 
stress axis; is the tangent modulus of the aforementioned asymptotic line. 
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The strain-rate sensitivity function k(0) was determined from the 
experimental constant-strain-rate stress-strain curves. It was found that 

the following function fit the data nicely for annealed aluminum in the 

-4 3 -1 

strain-rate range of 10 ~ 10 s . 


k(0) - 1 - B log 
8 



( 21 ) 


e 

In equation (21), 0 is a material constant and 0_ is a reference strain-rate. 

S K 

The expression (21) was also shown in the same report to be applicable to 

-4 1 -1 

mild steel in the strain-rate range of 10 ~ 10 s . 

In Wu and Yip [11] this form of constitutive equation was also used to 
investigate the strain-rate history effect by means of low-high and high-low 
strain-rate change loading sequences. 


if 




gr 


9 


3. CREEP OF ALUMINUM 1100-0 


The constitutive equation (13) and the new definition of intrinsic 

tine given by Eq. (IS) are now applied to describe the creep phenomenon. 

To achieve this objective, a careful control of the entire creep test is 

in order. In a conventional creep testing, a weight pan loaded with weights 

is gradually released to avoid impact on the specimen. The creep strain is 

then the additional strain developed over the initial strain as time elapses. 

The initial strain referred to is developed initially by the dead weights. 

The initial strain history is usually not recorded in the creep test. 

Due to this nature in experimentation, any theoretical work aiming at the 

description of the creep curve is at best an approximation. * 

In the present representation it is assumed that a creep test consists 

of two stages. The first stage is assumed to be a tensile testing at con- 

stant plastic strain rate 8 . During this stage the intrinsic time increases 

o 

from zero to z q . The second stage is the creep stage where the stress is 
kept constant at o* and the intrinsic time increases from z q to z as creep 
strain developcs. Thus, equation (13) becomes 
z z 


o* - E 


+ E 


'o 

r 

'z 


6(z - z' ) 


6(z - z' ) 


_d_6 

dz 




-*') de 


dz 


7 dz' 


de 

dz 


> r 

7 dz' + E 

*9 


-a(z-z') dO 


dz 


7 dz* 


( 22 ) 


o o 

in which the fuuction p(z) is given by equation (12). 


In the case of long term creep, the effect of the initial strain history is 
negligible and the stress-time history may be represented by the Heaviside 
step function. In the case of short time creep, the initial strain history 
is found to have significant effect on the subsequent creep behavior. 
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It should be notsd that during the first stage of loading for 
aluminum, the following equation must be satisfied 

- - ln(l + 8 k 0 ) (23> 

Op 0 0 

where 0 Is a material constant that appeared in equation (17) and Is 

determined from the quasi-static stress-strain curve; k ■ k(9 ); and 8 is 

o o o 

the plaotic strain at the end of the flrat stage when z - z . 

Figure 1 shows a typical initial strain history recorded during a 

conventional creep test. Approximating this initial strain history by a 

straight line, it is seen that at z ■ z , dO/dtl + 6 . Since stress 

o t z m z o 

o 

depends on the strain-rate and its history, the actual stress at z q is not 

the same as that calculated using equation (18) by assuming a constant o . 

o 

Therefore, the stress must be corrected at z q to agree with reality. This 
correction of stress will be accounted for in the subsequent discussion 
concerning the second stage of the creep test. It is remarked, however, 
that if the creep test is well-controlled so that the -irst stage of test 

e 

is truly a constant 6 q test. The correction in stress mentioned above is 
then not necessary. 

Equation (22) is the constitutive creep equation which may be simplified 
by letting 


R(s) - Ej^ 



e -o(*-*') 


de 

dz' 


dz' 


Thus, aqua don (22) reduces to 
o* » I(z) + R(z) + o° 


I(z) - C x e” az 


(24) 


(23) 


where 


(26) 
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and 



_L 



(<*40 )* 

o 



It nay ba shown by dlffersntlating (24) that 


(27) 


dR 

da 



(28) 


which leads to 



by uaa of Eq. (25). The constants appeared in the above eq< '"’.Ion are given 

by 



(30) 


(33) 


During the creep process, the Intrinsic tine is defined by 


dt - k c ($) |de| (32) 

where k (0) is the strain-rate sensitivity function under creep condition, 
c 

The function k,(8) is different from k(6) in general, where the latter 
function was defined fn (15) under constant-strain-rate condition. Since 
the creep behavior depends on stress as a parameter, the function k c (0) may 
also depend on stress as a parameter. 
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Using (32), (26) and (17) Equation (25) reduces further to 

k (6) ■ o° e®*/(o* - R(z) - C- e az ) (33) 

c y x 

Given a function k (6) at the strain rate range of creep, equation (33) can 
c 

thus lead to the determination of a creep curve. 

The strain-rate sensitivity function k (6) will now be discussed. It 

c 

is assumed that the functional form of equation (21) is still valid for 
creep. Thus, 

ye> - 1 - 8 C la ^ (34) 

in which 6 denotes the creep rate and the constant $ c depends on stress as 
a parameter. The dependence of $ c on stress should be further investigated 
experimentally. 

The creep curves are determined by use of the Newton method of finite 
difference. The algorithm used in the calculation is presented in Appendix 
(A). 

To compare the theory with experiment, two sets of experimental results 

are needed. They are the set of constant-strain-rate stress-strain curves 

and the creep curves. These two sets of data should be obtained for the 

same material with the same heat treatment and tested at the same constant 

temperature. The only experimental data that are available to the authors 

are those for Aluminum 1100-0. Creep tests were conducted at 150® C (300°F) 

» 

by Wang and Onat [10] and the constant-strain-rate stress-strain curves 

were recorded by Senseny, Duffy and Hawley [12] at 25°C (298°K) and 250°C 

-4 -1 2 -1 

(523°K) for shear-strain rates of 2 x 10 s and 3 * .10 s . To determine 
the material constants, a set of experimental constant-strain-rate stress-strain 
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curves at 150*C were estimated by linear interpolation from the curves of 
Ref. [12]. 

In the creep tests by Wang and Onat specimens were prepared from two 

different rods for the tests referred. Thus, different creep curves were 

obtained for the two groups of specimens even with the same heat treatment 

and same constant stress. The material properties for the two groups of 

specimens are therefore considered to be slightly different. The following 

* 

two sets of material constant have been determined: 


Specimens Group I 

0 R - 1.30 x 10’ 5 s -1 
• _4 — I 

e - 1.20 x 10 s 

o 

o - 4.622 x 10 2 
B - 4.444 

E x - 6.596 x io 3 MPa (9.566 x 10 5 ps i) 

6 - 3.643 x i(f 2 
8 

E q - 6.895 x io 4 MPa(1.0 xlO 7 psi) 

0 ° - 13.790 MPa (2000 psi) 


Specimens Group II 
0 R - 1.40 x io" 5 s" 1 
e - 1.20 x 10' 4 s -1 

O 

o - 5.714 •< 10 2 
B - 5.494 

E x - 8.154 x io 3 MPa (1.183 x 10 6 psi) 

B « 3.339 x io' 2 
s 

E - 6.895 x io 4 MPa (1.0 x io 7 psi) 

0 ° - 13.790 MPa (2000 psi) 


The theoretical and experimental results are shown in Figures 2-4 for 

specimens Group I and in Figures 5-7 for Group II. The values for fi and 

c 

the stress correction mentioned earlier are listed in Table 1 for each stress 
level. The results show that good agreement has been obtained for Group I 
specimens. For Group II, only the three cases with Intermediate stress 


The constants a, B, E, E^, Oy, B s and 0_ are determined from the constant— 
strain-rate stress-strain curves. The strain rate $ during the first stage 
of the creep test must be chosen so that 6 Q and z Q are consistent with 
equation (23). Thus, 0 O is the only parameter which needs to be assigned 
in the description of a conventional creep test. 
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levels have achieved reasonably good agreement using the constants listed. 
The agreement is not good for the two lover stress curves at small time and 
also for the curve with high creep stress. It should be noted that the 
present theoretical results have been obtained by assuming that only two 
Internal state variables are significant. To predict the creep curves at 
widely varied stress levels, an additional internal variable may have to be 
used. Whether this is indeed the case awaits further investigation. 

It may be concluded that for Aluminum 1100-0, creep depends on the 
effects of strain-hardening and strain-rate history. If 8^ is dependent on 
stress as mentioned above, the proposed theory can reasonably describe the 
creep curves with moderately varied stress levels. 
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4. CREEP RECOVERY 


In this section, creep recovery is discussed by means of the constitu- 
tive equation (13) with the Kernel function given by (12). The specimen is 
first subject to creep at stress o* for a duration of time and then the 
stress is reduced to a**. The integral in equation (13) is decomposed into 
four integrals each representing a stage of the creep-creep recovery phenom- 
enon. Thus, 


a** 


Z Z. Z, 2 

r-r-r-i 

*9 *9 *9 


(35) 


The first two integrals on the right-hand side of (35) describe the creep 
stage which was discussed in the previous section. The stress is then 
reduced from o* to a** when the intrinsic time reaches z^. The third inte- 
gral represents the instantaneous elastic response during unloading. Since 
the intrinsic time increases only with the plastic strain the value of this 
integral is zero. Finally, the fourth integral represents the recovery 
stage. 

The first integral may be integrated to give I(z). The second integral 
may be denoted by 


R x (z) 


- E 


r 


p(z - 


z’)ffrdz’ 


It can be found by differentiation that 


(36) 


dR r (z) 


dz 


-ort^z) 


(37) 
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whose soluClon is 



Rj^(s) ■ C e aZ 

(38) 

vaere 





R 1 (z 1 ) - R(zp 

(40) 

The fourth Integral is expressed as 



e <, f *<* - *•> w iz ' * If + v*> < 4i > 

z i 

where 

R 2 (z) - Ej f r a<s - *'> ffrdz' (42) 

*1 

Therefore, the constitutive equation for creep recovery may be rewritten as 
o** * I(z) + R^(z) + °y ff + R 2^ z ) (43) 

or 

~ - [o** - I(z) - R 1 (z) - R 2 (z)] (44) 

o 

y 

The function R 2 (z) may be found in closed form, which may be achieved 
by differentiating equation (42) to obtain 
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dKAz) .. 

~ — ’ E 1 ii - “*2 (,) 


( 45 ) 


Upon substituting (44) Into (45) , the resulting equation can be Integrated . 
The result Is 



where 


and 


C - C. + C 

X 1 



R 2 (*j) “ 0 


(47) 

(48) 


(49) 


Therefore, It Is obvious from equation (44) that ^ is a function of z 

£ felt 

only if a and o are specified. Using now equations (32) and (17), 
equation (44) can be written as 


** . 

k (8) 

c 


o Bz 
_X 


o** - I(z) - R. (z) - R,(z) 


(50) 


in which k denotes the function k (0) during the stage of recovery. The 
c c 

significance of sign "±" will be discussed later in this section. It is 

** 

noted from equation (50) that Is a function of z. 

A procedure similar to that used for creep stage is now followed. 
Assuming that the relation (32) holds also during the stage of recovery, 
the strain-rate at recovery may be expressed by 
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. . (l-k**(*))/8 

6(s) - 0 R e c (51) 

An algorithm of numerical calculation for creep recovery la presented in 
Appendix (B). 

The experimental results of Wang and Onat for Aluminum 1100-0 are again 
used for comparison. Three examples were considered in this calculation. 

The first (Test I) involved a specimen crept at constant stress of 24.13 MPa 
(3500 psl) for 4 hours and then unloaded to a stress of 1.38 MPa (200 psl) 
which was far below the yield stress of 13.79 MPa (2000 psl). In this cal- 
culation, the equation (50) took the minus sign since the normal strain 

** 

recovery effect was predominant at such a low stress of a ■ 1.38 MPa. 

The normal strain recovery effect will be elaborated later in this section. 

For this case, the material constants were determined by the method 

described in the previous sections as follows: 

0 - 1.20 x io“ 4 s" 1 
o 

e R - 1.40 x io ' 5 s" 1 
0 - 5.259 
a - 4.944 x 10 2 

E x » 7.06 x 10 3 MPa (1.024 x 10 6 psi) 

0 - 3.858 x 10" 2 
c 

Both theory and experiment are shown in Fig. 8 for comparison. 

The second and the third examples (Test II and III) involve specimens 
loaded at constant stress of 27.58 MPa (4000 psi) for one and four hours, 
respectively, and then unloaded to a stress of 25.86 MPa (3750 psi). In 
both cases, equation (50) took the positive sign since the forward creep 
effect was predominant at the stress specified (25.86 MPa). The details 




. t 


1 t 

t 
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concerning the forward creep and the recovery effects will be further 
discussed . 

For these two cases, the material constant were determined as: 

8 - 1.20 x 10’ 4 s" 1 
o 

e R ■ 1.40 x io 5 s” 1 
6 - 5.494 
o - 5.714 x 10 2 

E x - 8.16 x 10 3 MPa (1.183 * 10 6 psi) 

8 - 3.339 x i(f 2 
c 

Figure 9 shows the theoretical and the experimental results. It Is seen 
that the agreement Is quite good for a short time. The results clso show 
that the amount of forward creep during the process of recovery is signifi- 
cantly affected by the length of creep time before unloading. For the same 
creep condition, the earlier the unloading takes place the greater the for- 
ward creep rate. 

In Test II, the theoretical curve rises above the experimental curve 
four hours after unloading has taken place. Its Implication Is that the 
present theory using two internal state variables is not adequate in de- 
scribing creep recovery at larger time. For better accuracy additional 
terms (corresponding to more internal state variables) may have to be 
introduced in equation (12) . 

The nature of creep recovery will now be further discussed. Experiments 

show that a critical stress o cr exists. If during creep recovery, o > a cr 

then forward creep will occur, i.e., the recovery curve has a positive slope. 
** 

If o < o cr , then normal strofn recovery is found which corresponds to the 
curve with negative creep rate. This experimental finding has been reported 



i 


{ 

1 

j 

i 

i 

i 

4 
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In the literature (see, for instance. Ref. [1], p. 24), In particular the 

experiments or Wang and Onat [10] also confirmed thie observation. In the 

three examples mentioned earlier in this section, Test I is the case of 

normal strain recovery with o < o _ and Tests II and III are the caees of 

cr 

ft* 

forward creep with o > o • 

cr 

The observed phenomena of forward creep and normal strain recovery 

will now be explained using the present theory. Just before unloading, 

* - 
the stress is a and the intrinsic time is denoted by z^. From equation (25) , 

the following equation is obtained 


* 


o 


I(z") + R(z‘) + 0° 



(52) 


ft* + 

Immediately after unloading, the stress is o , the intrinsic time is z^, 

and from (43) the equation below is obtained 


** 


i(‘t> + h<$ * °y f 


(53) 


But I(z), R(zK R^(z), and R£(z) are all continuous so that at z^, 
R^(z*) • R(Zj) and I(z*) ■ I(Zj). Therefore, equations (52) and (53) 
can be combined to yield 


* ** o d0 

0 -° "°y Tz 


- a 


o d0 
y dz 


(54) 


Furthermore, by using equations (32) and (17), it may be found that 

o 


de 

dz 


S e pz i 

k c 


( 55 ) 
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«nd 


d8 

dr + 
*1 



(56) 


The "±" sign In (56) corresponds to the absolute value sign In equation (32) . 
In the case of forward creep, where 0 > 0, the positive sign should be used; 
whereas the negative sign is used when the case of normal strain recovery 
is being considered, i.e. , when 0 < 0. 

Substituting equations (55) and (56) into (54) , the following expres- 
sion for the strain-rate sensitivity function is obtained 



(57) 


* ** 

In the special case of no unloading, i.e., o - o , equation (57) leads 
a* 

to k ■ k (only the positive sign is meaningful due to positive creep 
0 c 

rate), which of course represents the case of creep. 

When at z^, equation (50) is reduced to 

Bz, 


k c <«!> • * 


o '.' 1 

-1 


“ ** , v 

o - I(z,) - R 1 (* 1 ) 


(58) 


Utilizing this equation, the following observations may be made. On the 

** 

right-hand side of equation (58), the denominator is zero when a - I(z^) 

+ R. (z.). Therefore, a critical stress o may be defined so that when 
lx cr 

a ■ o cr , k c (z^) approaches to infinity, which in turn gives 0 « 0 by use 

of equation (51). Thus, creep recovery at the critical stress a will 

cr 

result in a taro strain recovery rate at z^. For z > z^, the rate of 
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recovery le governed by equeeion (50) in which i 0. For ° > 

the denominator le greater than zero, thus the plus sign must be chosen in 

(58) for k** to be positive. As mentioned earlier, the plus sign corre- 

, ** 

•ponds to the case of forward creep. On the contrary, if a < a^, then 
the minus sign should be used, which of course corresponds to the case of 
normal strain recovery. 

A further observation related to equation (58) is that in the case of 

** ** 

forward creep, the greater rhe stress a is, the smaller is k £ ; hence, 
the greater is the rate of forward creep. In the case of normal strain 
recovery, the reverse is true, i.e., the smaller the stress a is, the 

HtUt 

smaller la k ; thus, the greater is the rate of recovery. The above 
c 

observation agrees with Garofalo [1] and Wang and Onat [10]. 
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5. STRESS RELAXATION 

A preliminary investigation of stress relaxation is reported in this 
section. Due to the non-availability to the authors of the stress relaxa- 
tion data for Aluminum 1100-0 at 150 # C (300 # F), only the theoretical results 
will be discussed. 

To treat the subject of stress relaxation using the present theory, it 
is recalled that 

de - de - I 2 (59) 

o 

During the process of stress relaxation, the total strain is kept constant. 
Thus, 

# 

e - - £- (60) 

o 

and the plastic strain-rate is seen to be related to the rate of change of 
stress during stress relaxation. 

The constitutive equation governing stress relaxation is again reduced 
from the general constitutive equation given by (13). The relaxation test 
is divided into two stages as in creep test. The first stage is that of 
the constant-strain-rate tension test and the second stage is the stage of 
stress relaxation. 

It may be shown that the following governing equation can be reduced 
directly from equation (13): 

o(z) - I(z) + R(z) + o° (61) 

where I(z) and R(z) were previously defined in equations (26) and (24). 


' * I 

i 
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The symbol z Is Again used to denote the Intrinsic time when the second 


o 


stage of test begins. 

For stress relaxation, the following equations are again obtained 


. E IS . * 

dz 1 d* 


(62) 


ff - J c ; [o(z) - I(z) - R(z) ] 
o 

y 


(63) 


These equations may be combined to yield a first order differential equation 

in R(z). The following solution for R is obtained by requiring that R(z ) -0: 

o 


v °y J 


e Az g(z‘) dz’> (64) 


where A and are defined in (30) and (27) and 


g(z) - -J o(z) 

°y 


(65) 


During the relaxation process, the intrinsic time Is defined by 


dc - k r (e)|de| 


( 66 ) 


where k f (6) is the strain-rate sensitivity function under the condition of 

relaxation. This function must be different from k (6) which is the strain- 

c 

rate sensitivity function for creep. During relaxation the total strain Is 

kept constant and this constraint should be reflected through k f (0) . It 

Is anticipated that k (8) and k (8) are related. However, no explicit 

c r 

relationship between these two functions can be obtained without further 
experimentation . 


1 

* 
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It is assumed that tha function (0) takas the following fora 



(67) 


where $ r may depend on the magnitude of plastic strain when relaxation stage 
begins. The fact that the theoretical results to be presented later agree 
qualitatively with the usual atress relaxation curves Indicates that the 
functional fora given by (67) Is reasonable. 

It may be obtained from (63) and (66) that 

0 ° . 6 ‘ 

k r (,> ‘ o(z) - I(t) - R(z) <68) 

In addition, equation (60) leads to 


o(t) • -^,6 r 


U-k r >/8 r 


(69) 


Thus, equations (68) and (69) govern the relaxation behavior of a material 
An algorithm for the numerical calculation of the relaxation curves are 
given in Appendix (C). The theoretical results are shown in Fig. 10. 
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6. CONCLUSIONS 

The endochronlc theory of viscoplasticity has been applied to discuss 
creep, creep recovery, and stress relaxation at the small strain and short 
time range. The following conclusions may be drawn from this investigation: 

(1) The governing constitutive equations for constant-strain-rate 
stress-strain behavior, creep, creep recovery, and stress relaxation have 
been derived by Imposing appropriate constraints on the general constitutive 
equation (13). 

(2) A set of material constants has been found which correlate strain- 
hardening, creep, creep recovery, and stress relaxation. 

(3) The theory predicts with reasonable accuracy the creep and creep 
recovery behaviors at short time. 

(4) The strain-rate history at prestraining stage affects the subse- 
quent creep. 

(5) A critical stress o cr has been established for creep recovery. 

itik 

If ^ > o , forward creep will occur. If o < o , then normal strain 

cr r cr 

recovery will take place. 

(6) The correlation between the st.ain-rate sensitivity functions 

k<0), k (6), and k (0) cannot be obtained without further experimental 
c r 

investigation. 
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APPENDIX (A) - Algorithm for the Computation 
of Creep Curves 


Based on equations (33) and (34), the creep curves are computed by use 

of the Newton method of finite difference. The following steps are suggested 

for computing creep curves obtained by the conventional method of creep 
* 

testing: 


Step 1: Determine the material constants a, B» E , E. , a , 8 , and 0_ from 

o x y s r 

the experimental constant-strain-rate stress-strain curves. 

* . * 

Step 2: Compute the plastic strain 6 ■ c - a /E where e and o are 

o o o o 

known values. 

Step 3: Try and error procedure is used to determine 0 q . For each 0 q chosen, 

C and z can be computed by use of equations (21), (32), and (23). 
o o 

The rest of the procedures are then followed to compute the creep 
curve. However, if the computed creep curve does not agree with 


the experimental curve, a new value of 6 q is assumed. 

Step 4: For each step, compute k^ n ^(0^) by substituting z^ and o* into 

c 

equation (33), where n starts from zero and is a positive integer. 

;(n) 


Step 5: Use equation (34) to compute 0 

Step 6: Choose a reasonably small step size At and compute A0 

.(n) 


(n) 


e (n) At 


and AC 

.( 0 ) 


k <n) Ae (n) 

C 


Step 7. Set t' 7 - 0, C 
„(n+l) 


( 0 ) 


V 2 


( 0 ) 


then compute 


t (n) + At 


* • g 

For well controlled creep test, i.e., 0 O • constant, the quantities 0 O and 
0 O are known. The intrinsic time z Q and the stress o* can be calculated 
from equations (21), (23), and (18). Noting that the total strain at the 
beginning of creep stage is e Q » 0 O + 0*/E o » the computation may proceed 
by following steps 4 through 8. 
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g(tt+l) m 9 (n) + ^gC®) 

e (n+l) . g(n+l) + £ (n) 
z (r&l) m c (n) + Ac (n) 


, (•*■!) - * (n) + A in 


ec (n) J 


The vs. t^ profile is the computed creep curve. 

Step 8: Specify a terminal time T. Ii is greater than T, stop 

the computation. Otherwise, return to Step 4 with t^ n \ 
e (n) , C (n) , and z (n) replaced by t (n+1) , 8 (n+1) , e (n+1) , C (n+1) 
and z (n+1) , respectively. 
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APPENDIX (B) - Algorithm for the Computation 
of Creep Recovery Curves 


Curves of creep recovery may be computed from equations (SO) and (51) 
using the following procedures: 


Step 1: 
Step 2: 
Step 3: 

Step A: 

Step 5: 
Step 6: 

Step 7: 
Step 8: 
Step 9: 


Use all previously determined material constants. 
Calculate the critical stress o 


cr 


I( Zl ) + R 1 (z 1 ). 


** ** 

Compare o with o . If o > ° cr » use positive sign in equation 

(50). Otherwise, the negative sign should be used. 

**(n) (n) ** 

For each step, compute by substituting z and a into 

equation (50), where n starts from zero and is a positive integer. 

• (n) 

Use equation (51) to compute 0 . 

Choose a reasonably small step size At and compute A0^ * 0^At 


(n) 


and AC 
Set t* 0 ^ 


c 




( 0 ) 


V 


.( 0 ) 


* z^ where t^ is the time at the 


( 0 ) * **si ( 0 ) 

onset of recovery. Compute e * - (o - o )/E o , where e 


and are the recovery and creep strain at t^, respectively. 

using the 
(n) (n) 


Compute t (n+1) , 0 (n+1) , e (n+1) , C (n+1 \ and z (n+1) 


equations in Step 7 of Appendix (A) . The e 
the computed creep recovery curve. 

Same as Step 8 of Appendix (A) . 


vs. t profile is 
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APPENDIX (C) - Algorithm for the Computation 
of Stress Relaxation Curves 


Equations (68) and (69) are used for the compuation of the stress 
relaxation curves. In the compuation of R(z) from (64), the function g(z) 
is assumed piece-wise linear. The following steps are followed: 


Step 1: Use same material constants as those in Appendix (A). 

Step 2: Compute 0 q , and z q using Steps 2 and 3 of Appendix (A). 

Step 3: When n * 0, set o^ » o r , 4^ * 4 , z^ * z , - 0, where 

o o o 

x 

n starts from zero and is a positive integer, o q is the initial 
relaxation stress. 

Step 4: For each step, except n * 0, compute R(z^) with the assumption 

of piece-wise linear o(z): 


o (n) - . (n-1). , 0-1) 

—7 (2 - z ) + a 

Az 


If n ■ 0, set R(z^ n ') * 0. 

Step 5: Compute by substituting z^ n \ and R(z^) into 

equation (68) . n starts from zero in this computation. 

• ( jx) 

Step 6: Calculate stress rate o from equation (69). 

Step 7: Choose a reasonably small step size At and compute the following: 

4o (n) . ; <n, 4t 

. , k (n> ; (n) 

._(n) _ . r 


Ac v,,/ 

1 + 64 
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Step 8: Compute 

t (n+l) . t („) + 4t („) 

0 (n «) . 0 (n) + M (n) 

t (n+1) . t <”> + at (n) 

2 (n+1) . z (n) + Az ( "> 

The a ^ ^ vs. t^ ^ profile is the computed stress relaxation curve. 
Step 9: Specify a terminal time T. If t^ n+1 ^ is greater than T, stop the 

computation, jcherwise, return to Step 4 with t^ n \ a^ n \ ^ n \ 
and z<"> replaced by o <n+1 » < ( " +1 \ s<” +1 >, respectively. 




TABLE 1 


Corrections in Creep for Various Creep Stresses 


Specimen 

Group 

Creep Stress 
MPa (psi) 

8 c 

Stress Correction 
MPa (psi) 

I 

24.13 

(3500) 

6 s 

-1.990 

(-288.68) 


25.86 

(3750) 

o.7i e s 

-1.643 

(-238.34) 

II 

17.24 

(2500) 


+0.454 

(+65.80) 


20.68 

(3000) 

8 s 

-0.869 

(-125.97) 


24.13 

(3500) 

6 s 

-1.724 

(-250.06) 


25.86 

(3750) 

0.79 S 8 

-1.628 

(-236.06) 


27.58 

(4000) 

0.43 B 

8 

-1.467 

(-212.70) 
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1. Initial Strain Hlatory in Creep Tost 

2. Cons tan t-Strain-Rate Stress-Strain Curves of Aluminum 1100-0 at 150°C 
(300°F) - Specimens Group I 

3. Strain-Rate Sensitivity Function k for Group I 

A. Creep Curves for Aluminum 1100-0 at 150°C (330 # F) - Specimens Group I 

5. Constant-Strain-Rate Stress-Strain Curves of Aluminum 1100-0 at 150°C 
(300°F) - Specimens Group II 

6. Strain-Rate Sensitivity Function k for Group II 

7. Creep Curves for Aluminum 1100-0 at 150°C (300°F) - Specimens Group II 

8. Creep Recovery Curve for Aluminum 1100-0 at 150°C (300°F) - Test I 

9. Creep Recovery Curves for Aluminum 1100-0 at 150 # C (300°F) - Tests II 
and III 

10. Theoretical Stress Re le t t e ws Curves 
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